Dipolar gases in coupled one-dimensional lattices 



Marianne Baucr^ and Meera M. Parish^' ^ 

^Cavendish Laboratory, J J Thomson Avenue, Cambridge, CBS OHE, United Kingdom 
^London Centre for Nanotechnology, Gordon Street, London, WCIH OAH, United Kingdom 

(Dated: February 21, 2012) 

We consider dipolar bosons in two tubes of one-dimensional lattices, where the dipoles are aligned 
to be maximally repulsive and the particle filling fraction is the same in each tube. In the classical 
limit of zero inter-site hopping, the particles arrange themselves into an ordered crystal for any 
rational filling fraction, forming a complete devil's staircase like in the single tube case. Turning on 
hopping within each tube then gives rise to a competition between the crystalline Mott phases and 
a liquid of defects or solitons. However, for the two-tube case, we find that solitons from different 
tubes can bind into pairs for certain topologies of the filling fraction. This provides an intriguing 
example of pairing that is purely driven by correlations close to a Mott insulator. 
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Unconventional superconductors such as the cuprates 
have stimulated much research on exotic pairing phenom- 
ena in low dimensions Superconductivity in these 
materials is typically observed upon doping a half-filled 
Mott insulating state, and thus there has been particu- 
lar interest in whether pairing can be purely driven by 
strong electron correlations near a Mott insulator. Model 
systems that have been used to investigate this effect are 
the spin ladders, where one-dimensional (ID) chains of 
electrons are coupled via electron hopping and magnetic 
exchange interactions 0, Here, when the ladder is 
doped, one can obtain pairing between holes on different 
chains depending on the number of chains or 'legs' in the 
ladder 0]. However, the focus thus far has been on Mott 
insulators derived from short-range, on-site interactions 
and so an intriguing and physically relevant question is 
how long-range interactions will affect the physics. 

It is well known that long-range interactions in ID can 
lead to exceptionally intricate crystalline ground states. 
This is most evident in the case of classical particles on 
a lattice interacting via repulsive infinite-range convex 
potentials 0, Q . Here, in the absence of any quantum 
kinetic energy, the particles arrange themselves into an 
ordered crystal, commensurate with the underlying lat- 
tice, for any rational filling fraction. Indeed, it can be 
shown that the filling fraction as a function of the particle 
chemical potential /i forms a complete devil's staircase, 
where every rational filling fraction between and 1 en- 
joys a region of stability within a finite interval of /x 
Perturbing away from the classical limit, one finds that 
kinetic energy destroys the complete devil's staircase, but 
signatures of the staircase structure remain in the form of 
Mott lobes 01 • This structure has also been predicted in 
recent studies of the opposite weak-coupling limit, where 
the lattice potential is weak Q. 

In this Letter, being motivated by the ladder systems, 
we extend the above problem to investigate the case of 
two ID lattices that are purely coupled by the long-range 
interactions. Such coupled ID systems have recently gen- 
erated much interest in the context of ultracold atomic 
gases owing to the possibility of confining polar molecules 
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FIG. 1: (Color online) Arrangement of dipolar bosons (filled 
circles) in two ID lattices separated by distance d, where 
the dipoles are directed out of the page, (a) Commensu- 
rate ground state for filling fraction p/q = 3/8 in each tube. 
The ground state pattern can be determined by collapsing the 
two tubes onto one tube (b) and then using Hubbard's algo- 
rithm [3 for a single tube. The distances Xi, yi correspond to 
the nearest and next-nearest neighbor distances, respectively, 
on the collapsed tube. 



with long-range dipole-dipole interactions in reduced ge- 
ometries @-ll|. Thus far, theoretical studies of tubes 
coupled by dipolar interactions have concentrated on the 
continuum limit, whe re ther e is no strong lattice poten- 
tial within each tube 
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opposite, strong-lattice limit. 

We consider dipolar bosons in two identical tubes of 
ID lattices, where the boson filling fraction is the same 
in each tube and the dipoles are aligned so as to be 
maximally repulsive (see Fig. [1} . We focus on hard-core 
bosons, but our results will also apply to dipolar fermions 
as we note below. In the absence of hopping between lat- 
tice sites, we find a complete devil's staircase like in the 
single-tube case, but where the width of the steps can 
now depend on the intertube separation. For a given fill- 
ing fraction p/q in each tube, we find that the character 
of the soliton excitations (or domain walls) in the com- 
mensurate Mott phase is governed by whether q is even 
or odd. Thus, the quantum melting of the Mott phases 
exhibits an unusual "odd-even" dependence as the inter- 
tube distance is varied. Crucially, when q is even, we find 
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(a) Even q 



(b) Odd q 



configuration satisfies \xi — Xj \ < 1 for all pairs i, j. The 
same applies to 2nd neighbor distances yi = xi + X2, 
y2 ~ X2 + X3 etc., where we have \yi — yj\ < 1, and so on 
for higher order fc-th neighbor distances. This yields, for 
example, the configuration in Fig.[TJb) for filling fraction 
3/4. For Hubbard's algorithm to apply, it is sufficient 
to assume that the interaction potential is convex, i.e. 
V{r + 1) + V{r — 1) > 2V{r), where r can correspond to 
any fc-th neighbor distance. 

To determine the ground state for the two-tube system, 
we first take d — so that all the particles are effectively 
collapsed onto one tube. We then obtain the CGS for a 
single tube with filling fraction twice that of each tube, 
i.e. 2p/q. When the two tubes are drawn apart slightly 
with d ^ 1, the particle configuration is unperturbed 
(the CGS has a finite energy gap) and the repulsion is 
clearly minimized if we assign every second particle to 
each tube (see Fig. [T]). This means that the configura- 
tion in each tube corresponds to the single-tube CGS for 
filling p/q, since each tube only contains the sets of 2nd, 
4th, 2fc-th neighbor distances of the collapsed tube, 
all of which satisfy the CGS constraint discussed above. 
Increasing d further, Hubbard's argument ensures that 
this configuration in the two tubes remains the ground 

state as long as the intertube potential V (^^Jr^~+'dPj 

The index a = {1,2} denotes the two different tubes, remains c onvex w ith respect to r^ . From an inspection 



FIG. 2: (Color online) The configuration of dipolar bosons 
in the hole-like soliton states for fillings p/q = 1/2 (a) and 
p/q =1/3 (b), which have even and odd q, respectively. Top: 
configuration around one defect. Bottom: the arrangements 
of the defects along the tubes. 



that solitons on each tube bind together in pairs, an ef- 
fect that is driven by the topology of the crystalline phase 
rather than by an attractive interaction like in Ref. [l3| . 

We consider the following extended Bosc Hubbard 
Hamiltonian for two tubes: 



H = (cI+l,aCj,a +h.C.j + ^ V{rij) fli^ahj^a 

i,a. i>j,a 



fii,a = c|c(Ci.Q gives the density at site i on tube a, t is 
the hopping between neighboring sites, and V{r) =Vo/r^ 
is the dipole-dipole interaction with Vq proportional to 
the square of the dipole moment. The second and third 
terms in Eq. ((1} correspond to the intratube and in- 
tertube interactions, respectively. The interparticle dis- 
tances arc parameterised in terms of the horizontal dis- 
tance Vij ~ I * — j I , where we have set the lattice spacing 
in each tube to be 1, and the intertube distance d (see 
Fig. [1]). The on-site repulsion is taken to be infinite, 
so that we have hard-core bosons (or, equivalently, spin- 
less fermions 15|)- Note that the chemical potential [i is 
the same for each tube, corresponding to an equal filling 
fraction p/q in each tube. We consider filling fractions 
p/q < i so that the total energy of the system remains 
finite in the limit d — 0. However, we expect our main 
conclusions to also hold for ^ < p/q < 1. 

We begin by considering the classical case of zero hop- 
ping (i = 0). In the limit d — > 00. we clearly have two 
isolated tubes and thus we recover the classical commen- 
surate ground state (CGS) for a single tube [1,[H. Here, 
we have a crystal commensurate with the lattice: for 
density p/q in each tube, the configuration of particles 
has period q with p particles arranged in each period 
(assuming that p and q have no common factors). The 
configuration for a single tube can be determined using 
an algorithm proposed by Hubbard [3| • If we denote the 
distance between particles 1 and 2 as xi, the distance 
between particles 2 and 3 as X2 and so on, then the CGS 



of d^V (V?"^ + d^), we expect that convexity could be 
violated when d > 2ry . This could lead to a scenario 
where the crystal patterns in each tube become shifted 
with respect to one another once d > min(2rij). How- 
ever, we find that this shift does not occur for infinitely 
long tubes because of the terms in the energy involving 
large fc-th neighbor distances where convexity is not vi- 
olated. Thus, the CGS configuration remains the same 
for arbitrary d. Note, though, that this convexity viola- 
tion does have consequences for the defect state discussed 
below. 

To assess the stability of the two-tube CGS for each 
filling, we must consider adding/removing one particle 
to/from each tube. For the single-tube case, such an 
addition (removal) leads to the formation of solitons or 
defects in the tube pattern corresponding to an arrange- 
ment where the length of the unit cell is reduced (in- 
creased) by one vacant site 0, @|- Indeed, for filling frac- 
tion p/q, the added particle or hole fractionalizes into q 
solitons that are particle- or hole-like, respectively. These 
g-soliton states (gSS) can also be deduced from Hub- 
bard's algorithm 

Returning to the two-tube case, we determine the qSS 
by starting from the collapsed tube (d = 0) once again. 
By removing two particles (one from each tube) from the 
CGS for filling fraction 2p/q, we immediately obtain two 
different scenarios: for odd q, we generate 2x q solitons, 
while for even q, we have 2 x q/2 = q solitons since the 
crystal period in this case is q/2. Now, when we separate 
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the tubes, keeping d ^ 1, we obtain the single-tube hole- 
qSS for fining fraction p/ q on each tube, which gives q 
solitons on each tube. However, the difference in topol- 
ogy between odd and even q means that for even q, these 
solitonic defects are bound in pairs across the tubes, while 
for odd g, all the solitons repel each other and are max- 
imally spread out — see Fig. [2ja) and (b), respectively. 
The same situation applies to the particle-like solitons 
generated when a particle is added to each tube. 

The pairs of solitons for even q are "maximally bound" 
as in Fig. [5Ja) as long as the configuration correspond- 
ing to the qSS for the collapsed tube remains the lowest 
energy state. However, as discussed earlier, this is not 
guaranteed once d > min(2rij) and indeed we find that 
soliton pairs start to separate above a critical d, becom- 
ing more weakly bound as d is increased (see Fig. [3]). The 
solitons eventually unbind in the limit c? — >■ oo , where we 
recover the single-tube case. 

Following Ref. we determine the interval in chemi- 
cal potential A/u over which a given CGS is stable by com- 
paring the CGS energy with the energies of the particle- 
qSS and the holc-qSS. Like the single-tube case, we find 
that the two-tube system in the thermodynamic limit 
exhibits a complete devil's staircase, where every tube 
filling fraction p/q enjoys a region of stability A/i which 
all add up to fully cover the range of p.. However, we 
also find that the energy gap A/Lt depends on the inter- 
tube distance d when q is even, but not when q is odd, 
as follows: 

A//odd = ^ mqF{inq, 0) 

771 — 1 



h l-l 



(2m- l)q 



1=1 m=l 



-2y(Va2 + d2) + F(v/(a-l)2 + d2) + 



where F{a, d) 

V{^yia+ly + d^). Thus, for large q, Afi odd scales ap- 
proximately as l/f/'* for all d, while A/icvcn scales as 
(2/g)'* when d — > and as l/q"* when d -> oo. This 
is consistent with the odd-even behavior of the solitons 
in the two-tube system. 

We now turn to the case of finite hopping t > 0. Here, 
the solitonic excitations can melt the CGS Mott phase 
into a Luttinger liquid 01 , leading to a series of Mott 
lobes reminiscent of the phase diagram for the Bose Hub- 
bard model [3]. The phase boundary of each Mott lobe 
in the ii-t phase diagram (Fig. [J]) corresponds to when 
the energies of the CGS and its adjacent qSS are equal, 
i.e., for each t and d, we take EcasilJ') = Eqss{^J') |2a |. 
Similarly to Refs. we determine the energies using 

a strong-coupling expansion in t/Vo. We expand up to 
second order in t/Vo, which is sufficient to capture the 
basic shape of the Mott lobes for odd q like in Fig. SJ^a) , 




FIG. 3: (Color online) The separation h between bound 
particle-like solitons as a function of intertube distance d for 
p/q = 1/4 on each tube. The solitons move apart in steps of 
q = 4 and we define ft = to be the maximally-bound con- 
figuration (inset). The solitons start separating once d > 3.8. 



but is not always enough for even q as we discuss below. 
This pcrturbative approach is also never accurate at the 
lobe tip, where the energy gap of the CGS goes to zero, 
but we are in any case more concerned with the sides 
of the Mott lobe, where the CGS is doped with parti- 
cles (upper boundary) or holes (lower boundary). The 
strong-coupling expansion for the CGS energy possesses 
the standard form: the first-order correction is zero, while 
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the second order correction is E^Qg — —^u ^^^^ JSW] 
where N is the number of particles and AEi is the differ- 
ence in energies for i = between the ground state and 
the excited state created by hopping the z-th particle. 

The qSS is highly degenerate in the classical limit and 
thus the relevant states to consider in the perturbation 
theory are the momentum eigenstates of the soliton. For 
odd g, the solitons are unbound and thus the corrections 
to the qSS energy have a similar form to those for the 
single tube in Ref. 0]. In this case, the first-order cor- 
rection per added particle/hole is — 2ig, corresponding to 
solitons with zero momentum (which lie at the bottom 
of the band). Figure HJ^a) depicts the Mott lobes up to 
second order int/Vb ioi p/q = 1/3, representative of the 
behavior for odd g as a function of d. The width of the 
lobes at t = corresponds to A/^odd- We see that the 
shape of the lobe changes very little with d, owing to 
the fact that the solitons have the same structure in the 
limits d — > and d — >■ oo . 

For even q and d ^ 1, the first-order correction to 
the qSS energy is instead —tq per added particle/hole, 
since the solitons are now bound in pairs. Here, each 
maximally-bound pair is hopped by q/2 sites when a 
soliton in the pair is hopped by q sites (see Fig. [2][a)). 
The second-order correction for the qSS (with one parti- 
cle/hole added to each tube) has the form 
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N/2q 
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Like before, AEr-^a.p is the potential energy cost for hop- 
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FIG. 4: (Color online) Mott lobes in the n-t phase diagram corresponding to tube filling fractions (a) p/q =1/3 and (b) 
p/q = 1/4 for different values of intratube distance d. In the limit d ^ 0, the lobes for p/q = 1/4 and p/q — 1/3 become 
equivalent to the lobes for a single tube with filling fraction 1/2 and 2/3, respectively. The inset in (b) zooms into the Mott 
lobes for larger distances. 



ping a particle, where is the position of the hopped 
particle on tube a with respect to the soliton on its 
tube. The particle can be hopped towards (+) or away 
from (— ) its soliton. Note that degenerate states, where 

(2) 

Ai?ri,ct.^ = 0, arc excluded. The first term of E^gg is in- 
dependent of a and corresponds to the soliton pair prop- 
agating via an intermediate state where the soliton pair 
pulls apart. 

Beyond a critical d, the solitonic pairs at i = start 
to separate and the first-order correction to the qSS en- 
ergy can then be zero. This suggests that our simple 
perturbative expansion is inadequate for large d since we 
expect the first-order correction to be —2tq per added 
particle/hole in the single-tube limit d ^ oo. Indeed, 
we find that the energy difference between states with 
different soliton separations (Fig. [3]) rapidly approaches 
zero with increasing d so that the second-order correction 

(2) 

^qSS becomes sizeable even for relatively small d. This 
implies that we must include states \h) with different soli- 
ton separations h in our strong-coupling approximation 
for the energy. We construct the effective Hamiltonian 
j^eff ^Yje set of states {\h)} by including any contri- 
butions of excited states not included in {| ft,)} up to sec- 
ond order in t/Vo. For instance, the diagonal elements 
of will contain the potential energies for each \h) 
plus a t'^ term due to fluctuations of the underlying com- 
mensurate pattern. This approach effectively amounts to 
performing a Schricffcr- Wolff transformation [3, [l^ on 
the original Hamiltonian ([T]). We then diagonalise H^^ to 
determine the energy of the qSS. In practice, we consider 
{\h)} up to a maximum separation /imax, where /imax is 
large enough that the lowest eigenvalue of H"^ does not 
depend on /imax- 

FigureHI^b) depicts the resulting Mott lobes for even q, 
where we consider p/q = 1/4 for various d. For d > 2, we 
determine the lobes using H"^ from the Schricffcr- Wolff 
transformation. We see that the size of the Mott lobes 



changes dramatically with d and is dominated by how the 
width A/iovon at t = depends on d. Note, also, that the 
slope of the Mott boundaries near t = becomes steeper 
with increasing d, approaching d^/dt = ^2q for top and 
bottom boundaries, respectively. The evolution of the 
Mott lobes with d provides an unambiguous signature of 
the soliton pairing for even q. 

In principle, one can determine the nature of the 
Mott transition by approaching it from the Luttinger liq- 
uid [ll,!!^. For a single tube, the transition at the Mott 
lobe tip is of the Kosterlitz-Thouless type, while the tran- 
sition everywhere else is described by a two-band model 
of quasiparticles which are gapped in the Mott phase 0] ■ 
Further work is required to ascertain how this scenario 
is affected by soliton pairing and the different soliton- 
soliton interactions in the two-tube case. 

Our predicted phase diagram should be experimentally 
realizable with ultracold dipolar atoms or molecules. For 
a typical lattice spacing of 500nm within each tube and 
tunable lattice depths of 5Er to 30Er, where Er is the 
recoil energy, the hopping t can range from O.li?,- to 
5 X 10-'^ Er. Thus, for ''°K*^Rb polar molecules with 
dipole moment ^0.2 Debye as in current experiment ll| . 
we obtain t/Vo — 0-02 — 3, which is sufficient to observe 
the behavior of the Mott lobes in Fig.djb). To enhance 
Vq, one can consider Rydberg-dressed atoms, where the 
effective dipole moment can reach 10 Debye. This al- 
lows one to lower t/Vo to well within the Mott lobes of 
Fig. m The commensurate Mott phases could be probed 
by Bragg spectroscopy [2l| and, for the case of Rydberg- 
dressed atoms, the solitonic pairs could be detected by 
single-atom- resolved fluorescence imaging [2^ . 

We finally note that our two-tube system is also po- 
tentially connected to quantum Hall bilayers since the 
strong-coupHng Hmit of each ID lattice corresponds to 
the Tao-Thouless limit of the fractional quantum Hall 
effect [13. 
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